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Abstract. We quantify the uncertainty in weak lensing 
mass estimates of clusters of galaxies, caused by distant 
(uncorrelated) large scale structure along the line of sight. 
We find that the effect is fairly small for deep observa- 
tions (20 < R < 26) of massive clusters (a = 1000 km/s) 
at intermediate redshifts, where the bulk of the sources are 
at high redshifts compared to the cluster redshift. If the 
lensing signal is measured out to 1.5 Mpc the typical 
1(7 relative uncertainty in the mass is about 6%. How- 
ever, in other situations the induced uncertainty can be 
larger. For instance, in the case of nearby clusters, such as 
the Coma cluster, background structures introduce a con- 
siderable uncertainty in the mass, limiting the maximum 
achievable signal-to-noise ratio to ~ 7, even for deep ob- 
servations. The noise in the cluster mass estimate caused 
by the large scale structure increases with increasing aper- 
ture size, which will also complicate attempts to constrain 
cluster mass profiles at large distances from the cluster 
centre. However, the distant large scale structure studied 
here can be considered an additional (statistical) source 
of error, and by averaging the results of several clusters 
the noise is decreased. 

Key words: cosmology - gravitational lensing - galaxies: 
clusters: - dark matter 



1. Introduction 

Light rays are deflected as they travel through the uni- 
verse, and as a result our view of the distant, faint galaxies 
is somewhat distorted. Measuring these lensing induced 
distortions has proven to be a valuable tool in observa- 
tional cosmology (see the review by Mellier 1999): the lens- 
ing signal is sensitive to the projected mass surface density 
of the deflector, and through lensing one can probe the 
mass distribution without having to rely on assumptions 
about the dynamical state and nature of the deflecting 
matter. 
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When the surface density of the lens is sufficiently high, 
the distortions can be so large that multiple images can be 
seen, which provide strong constraints on the mass distri- 
bution in cluster cores. In this paper we will concentrate 
on the weak lensing regime, which allows us to study the 
cluster mass distribution on larger scales. 

Weak lensing is a powerful tool to probe the distribu- 
tion of matter in the universe on various scales. One ex- 
citing application is the measurement of the lensing signal 
induced by large scale structure (e.g. Blandford et al. 1991; 
Kaiser 1992; Bernardeau, van Waerbeke, & Mellier 1997; 
Schneider et al. 1998, SWJK hereafter), and recently sev- 
eral groups have reported detections of this signal (Bacon 
et al. 2000; Kaiser, Wilson, & Luppino 2000; van Waer- 
beke et al. 2000; Wittman et al. 2000). 

The first successful weak lensing results were obtained 
by studying the systematic distortion of the images of 
faint galaxies behind rich clusters of galaxies (e.g. Tyson, 
Wenk, & Valdes 1990). Measuring the masses and mass- 
to-light ratios of rich clusters is an important application 
of weak lensing, because comparison of the weak lensing 
result with other mass estimators (e.g., X-ray, S-Z effect, 
or velocities) can give important clues about the geome- 
try and dynamical state of the system (e.g, Zaroubi et al. 
1998). 

Weak lensing also has some drawbacks. First of all, 
the amplitude of the lensing signal depends on the redshift 
distribution of the faint background galaxies. These galax- 
ies are generally too faint to be included in spectroscopic 
redshift surveys, but photometric redshifts are thought to 
provide ample information on the redshift distribution of 
these galaxies (e.g., Hoekstra, Franx, & Kuijken 2000). 
Also the observed weak lensing signal is a measure of a 
surface density contrast, and one of the problems associ- 
ated with this is the mass sheet degeneracy (Gorenstein, 
Shapiro, & Falco 1988). One could try to get around this 
problem by measuring the lensing signal out to large dis- 
tances from the centre, and arguing that the cluster sur- 
face density should vanish there. 

Another potential problem is the fact that the lensing 
signal is sensitive to all matter along the line of sight. 
In general, weak lensing analyses of clusters assume that 
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the signal induced by the cluster under investigation is 
so dominant that the contribution from other structures 
along the line of sight can be neglected. 

The effect of structures, close to the cluster, such as 
the filaments that connect the cluster to the 'cosmic web', 
or an elongation along the line of sight, have been stud- 
ied through N-body simulations (e.g, Cen 1997, Metzler 
et al. 1999). Local structures introduce both a bias and 
an additional uncertainty in the cluster mass estimate, in 
particular when one for instance tries to measure the mass 
within r2oo- 

In this paper we study how distant structures affect 
weak lensing mass estimates of clusters of galaxies. Unlike 
the local structure, these are uncorrelated with the cluster, 
and they can be filaments, galaxy groups, or even other 
clusters. That such structures affect the mass estimates is 
easily seen in the case of a massive high redshift cluster, 
with a low redshift group in front of it. Assume the cluster 
is at a redshift of unity, with a group at say z = 0.3, and 
that the mass of the group is one tenth of that of the 
cluster: the contribution of the group to the total lensing 
signal is about one third, compared to the 10% in mass. 

In this paper we will quantify the uncertainty in weak 
lensing mass estimates of rich clusters caused by interven- 
ing matter. To do so we examine the predicted statistical 
properties of the contribution by large scale structure. We 
use an approach similar to what is done to predict the 
amplitude of the cosmic shear signal, which is used to 
constrain the cosmological parameters. Given a cosmolog- 
ical model, we estimate the uncertainty in cluster masses 
derived from weak lensing analyses. 

The structure of this paper is as follows. In section 2 we 
discuss two popular methods to measure cluster masses. 
In section 3 we derive the equations that describe the ef- 
fect of large scale structure on such mass estimates. The 
effect of large scale structure on mass estimates for clus- 
ters is quantified in section 4. We examine the dependence 
on cluster redshift, aperture size and source redshift dis- 
tribution. 

2. Estimating masses from weak lensing 

Cluster masses can be measured in several ways. In this 
paper we will use aperture mass statistics, which provide 
a filtered measurement of the of the surface mass density 
in an aperture (Schneider 1996, S96 hereafter): 



M ap (9)= [ e d 2 9'u(\e'\)K(e') : 

Jo 



(1) 



where U{9) is the weight or filter function. Provided that 
U{9) is compensated, i.e. 



db bU(b) = 0, 



(2) 



the aperture statistic can be related to the tangential shear 
7t- Unfortunately cannot be observed directly, but one 



can measure the distortion gx = 7t/(1 — k)- In the weak 
lensing regime, where k w 0, one can use jt ~ <7t, which 
is what we do throughout this paper. Then we have 



M ap (9)= f d 2 9Q{\9\) lT {9). 
Jo 

Here, Q{9) is related to U{9) through: 



Q(b) 



b 2 



db'b'U{b')-U{b). 



(3) 



(4) 



S96 has shown that the signal-to-noise ratio of the 
aperture mass measurement is optimal when <2(|0|) is the 
same as the (expected) tangential shear as a function of 
\9\. 

We want to quantify the uncertainty in mass measure- 
ments of clusters, and therefore require that the aperture 
statistic has a clear physical meaning. Studies of cosmic 
shear (see SWJK) use different choices for Q(|6>|). Below 
we discuss two statistics that have been used to estimate 
cluster masses. 



2.1. SIS model jit 

A popular model to describe the cluster mass distribution 
is the singular isothermal sphere (SIS) 



(5) 



where r E is the Einstein radius. The Einstein radius (in 
radians) is related to the velocity dispersion through 



r E = 4tt - (3 



(6) 



The tangential shear as a function of radius is simply given 

by 



7T (r) = n(r) = 



r_E_ 
2r ' 



The best fit SIS model is obtained by minimizing 



E 



7T,i - r E /(2n) 



(7) 



(8) 



with respect to r E , where jT,i are the individual measure- 
ments of the lensing signal at radius rj. The uncertainty 
in the measurement of the tangential shear, a 1 is given by 



o- 2 = 



7 27rnr i Ar i ' 



(9) 



with CTg a i the uncertainty in the measurement of an indi- 
vidual galaxy, n the number density of sources, and Ar^ 
the width of the ith radial bin. If we take the widths of all 
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the bins the same, and assume that the number density of 
sources is uniform, minimizing with respect to te gives 



(10) 



For a continuous profile the summation can be replaced 
by an integral. As mentioned before one measures the dis- 
tortion and not the shear, and equation [To] only works 
well when k is small. Also substructure in the cluster core 
changes the average tangential distortion (e.g., Hoekstra 
et al. 2000). To come around these problems one has to 
choose an inner radius for the fit, such that equation [l^ is 
a good approximation. 

If we assume that the tangential shear has been mea- 
sured from a distance i?o out to R, the Einstein radius 
satisfies 



te 



\n(R/R ) 



f 



dr Jt(t)- 



(11) 



Comparison with equation |3J shows that the best fit 
SIS model is an aperture mass measurement where Q{\0\) 
is given by 



Q{\0\) = 



(n9\n(R/R )) 




1 ,Ro< \0\ < R 
, elsewhere 



(12) 



The corresponding weight function for the surface den- 
sity U(8) is given by 



u(\e\) 



i 



n\n(R/R Q ) 



i/\e\- 



-2/R 
2/R 







o < \e\ < R 

R <\0\<R 
\9\ > R 



(13) 



The mass profiles of many clusters that have been stud- 
ied to date are consistent with a SIS model. Thus the 
weight function suggested above is close to the optimal 
choice. However, this choice for [/(I^Q is not a continuous 
function. A functional form which is continuous, yet simi- 
lar is suggested by S96, who discusses the use of aperture 
mass statistics to find mass concentrations in deep images. 
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Fig. 1. (a) The function gsig(l,0) as a function of W for 
a = 0.15, and = 200" (b) gf(l, 0) for a = f . Both curves 
show similar behaviour, with g^(l, 0) peaking at a similar 
values of 10. 



2.2. (statistic 



This gives 



Another statistic which has been widely used is the (- 
statistic (Fahlman et al. 1994): 



C(0i,d a ) = «(< 0x)-R{6i < \9\<6 2 ). 



(14) 



It gives the difference of the average dimensionless surface 
mass density in a circular aperture of radius 0\ and that in 
an annulus with 0\ < \0\ < 02- The corresponding weight 
function U(\0\) is given by: 



u(\e\) 



-{<ol-0\))- 



,\o\ < 0i 

,01 < \0\ < 02 

, elsewhere 



(15) 



Q(\o\) = 



0l{^(0l-0\))^ ,0 1 < \0\<02 







, elsewhere 



(16) 



Thus the (^-statistic can be related to the tangential shear 
through: 



Wu02) = -J% P d\n0( lT )(\9\). 

U 2 ~ U l JB X 



(17) 



The SIS model fit has only a simple physical meaning 
when the cluster mass profile is indeed isothermal. The 
interpretation of the ^-statistic is independent of the mass 



4 



H. Hoekstra: Effect of LSS on mass estimates 



distribution, although the statistic is not optimal in terms 
of signal-to-noise ratio. 

Under the assumption that one can estimate or mini- 
mize k{6\ < \9\ < 6 2 ), the average surface density in the 
control annulus, one obtains a measurement of k(< 9\), 
and thus of the cluster mass. One way to do so is to assume 
a model for the cluster surface mass density which is used 
to estimate the average surface density in the annulus (e.g. 
Hoekstra et al. 1998). This way the cluster mass estimate 
is less dependent on the assumed model compared to the 
SIS model fit. 



3. Contribution by LSS 

Under the assumption that the cluster is the only lensing 
structure in the field, the aperture mass M Q b s provides 
a direct measure of the cluster mass. However, the weak 
lensing signal is sensitive to all matter along the line of 
sight, and other structures, which we will refer to as large 
scale structure, will introduce an additional signal: 



M ohs {9) = M cluster (0) + M LSS (9). 



(18) 



The aperture mass statistic measures a density con- 
trast, and the expectation value for the contribution by 
uncorrelated large scale structure (A/lss) = (see also 
SWJK). Thus on average, distant structures do not intro- 
duce a bias in cluster mass estimates, but an uncertainty 
in the cluster mass of (-^lss) 1 ^ 2 - ^ s a resm t t ne effect 
of distant large scale structure can be treated as an addi- 
tional (statistical) source of error. Its relative importance 
can be decreased by averaging the results for several clus- 
ters, i.e. measuring their average mass. We note that on 
certain scales Mlss might be distributed skewly. 

The effect of 'local' (or correlated) large scale struc- 
ture has been studied using N-body simulations (e.g., Cen 
1997; Metzler et al. 1999). In this case either projection 
effects, or filaments connected to the cluster bias the weak 
lensing mass estimates towards somewhat higher values. 
This becomes particularly important when one tries to es- 
timate the enclosed mass within several Mpc's from the 
cluster centre (Cen 1997). 

SWJK have shown how one can relate the aperture 
mass statistic to the power spectrum of density fluctua- 
tions. Measuring {M 2 ^) 1 / 2 for a number of random fields 
allows one to constrain the projected power spectrum. In 
this paper the problem is reversed: given the power spec- 
trum, what is the resulting (M^g) 1 / 2 ? A detailed discus- 
sion about the use of aperture mass statistics and cosmic 
shear can be found in SWJK and Bartelmann & Schneider 
(1999). 

SWJK showed that the effect of large scale structure 
on the variance in the aperture mass statistic is given by: 



(M 2 ss )(9) = 2tt 



dl lP K (l)g(l, Of 



(19) 



where the effective projected power spectrum P K (l) is 
given by: 



BUS*. r dw (mXpJ-U; W )(20) 



4c 4 Jo 



a(w) 



Here w is the radial coordinate, a(w) the cosmic scale fac- 
tor, and fx(w) the comoving angular diameter distance. 
W(w) is the source averaged ratio of angular diameter dis- 
tances Di s /D s for a redshift distribution of sources p w (w): 



Jk{w') 



(21) 



The function <?(/, 9) in equation [T^ depends on the filter 
function f/(|0|), and is given by: 



9(1, 0) 



d4 cj>U(4>)J Q (i(t>). 



(22) 



As has been demonstrated in Jain & Seljak (1997) and 
SWJK it is crucial to use the non-linear power spectrum 
in equation 20. This power spectrum can be derived from 
the linearly evolved cosmological power spectrum follow- 
ing the prescriptions of Hamilton et al. (1991), and Pea- 
cock & Dodds (1996). 

3.1. Filter functions 

The form of the function g 2 (l, 9), which probes the pro- 
jected power spectrum, depends on the choice of U(\9\). 
For the SIS model fit we take Rq — aR, which yields 



gsisih 9) 



1.0 



2 {Jx(ald) ~ Ji(W)} + J dsJo(s) 

aW 



(23) 



We will use a typical value for a = 0.15 (e.g., Hoekstra et 
al. 1998). Larger values are used when much substructure 
is found in the cluster core (e.g., Hoekstra et al. 2000). 
The resulting ggigC) 0) is presented in figure [l]a. 
For the (-statistic we find 



J x (od6) -aJx{W) 
iraW(l - a 2 ) : 



(24) 



where we have taken the inner radius 9\ — a9. A typical 
value for a used in weak lensing mass estimates of clusters 
is |, i.e., the control annulus runs from the boundary of 
the region in which we want to determine the mass out to 
1.5 times this radius. Throughout this paper we will use 
this value for a. In figure |l]b the filter function g 2 (l, 9) is 
shown (taking a = |). 
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Fig. 2. Photometric redshift distributions from the HDF 
north and south. The counts have been normalized to a 
peak value of unity. The solid line is for galaxies with 
20 < R < 22, the dashed line is for 20 < R < 24, and 
the dotted line for 20 < R < 26. These distributions are 
used to compute the lensing signals and the contribution 
of large scale structure. 
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3.2. Photometric redshifts 

The projected power spectrum P K (k) depends on the red- 
shift distribution of the sources, and consequently we need 
to know, or assume, a redshift distribution for the sources. 
We will use the photometric redshift distributions in- 
ferred from the northern and southern Hubble Deep Fields 
(Fernandez-Soto, Lanzetta, Yahil 1999; Chen et al. 1998), 
which generally work well (e.g., Hoekstra, Franx, & Kui- 
jken 2000). There is, however, some field to field variation 
in the redshift distribution the sources, but we assume 
that the distributions obtained from the HDFs are rep- 
resentative for the whole sky. We note that most of the 
lensing by large scale structure is caused by structures at 
intermediate redshifts, and as a result the signal induced 
by the LSS is fairly robust against variations in the source 
redshift distribution. 

Using the colours of the galaxies in the HDFs we have 
computed their R magnitudes, which allows a direct com- 
parison with ground based observations. We use galax- 
ies fainter than R = 20. The resulting photometric red- 
shift distributions for different limiting magnitudes are 
presented in figure 



Fig. 3. (a) Scatter in the estimate of the Einstein radius 
Te due to large scale structure for the SCDM model (thick 
lines) and the OCDM model (thin lines). We use three 
different redshift distributions for each model. These are 
based on photometric redshifts inferred from both HDFs 
taking different limiting magnitudes. The solid line corre- 
sponds to 20 < R < 22, the dashed line to 20 < R < 24, 
and the dotted line to 20 < R < 26. We have indicated 
a radius of 1 h$Q Mpc for clusters at various redshifts. 
Panel (b) shows the results for the (^-statistic as a func- 
tion of aperture size. 



3.3. Aperture statistics 

Given the photometric redshift distribution discussed 
above we calculate the contribution of the LSS to the 
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n A 


C8 


SCDM 


0.5 


1.0 


0.0 


0.5 


OCDM 


0.7 


0.3 


0.0 


0.85 


ACDM 


0.7 


0.3 


0.7 


0.90 



Table 1. Parameters of the cosmological models consid- 
ered here. All three models are cluster normalized. 
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Fig. 4. (a) Relative importance of lensing by distant large scale structure on the measurement of the Einstein radius 
Te as a function of cluster redshift, when the shear is measured in an annulus from 0.225 to 1.5 Mpc. (b) the 
scatter in the measurement of te caused by large scale structure, (c) The observed Einstein radius as a function of 
cluster redshift. (d) same as panel (a), but now when the C — statistic is used to estimate k(< 1 h$Q Mpc). We have 
assumed that the lensing signal is measured out to 1.5 Mpc from the cluster centre, and that the average surface 
density R{1 < r < 1.5 h7 Q Mpc) is determined from a SIS model fit to the tangential shear, (e) scatter in R because 
of large scale structure, (f) k(< 1 Mpc) as a function of cluster redshift. The thick lines correspond to the SCDM 
model, and the thin lines to the open model. The values for either te and k turn out to be almost identical for these 
two cosmologies. 



measurements of te and the ^-statistic as a function of 
aperture size. We have considered three different cosmolo- 
gies, which are listed in table ^ 

Figure ||a shows the scatter in the value of the Einstein 
radius te as a function of aperture size for the SCDM 
(thick lines) and the OCDM (thin lines) models. For clar- 
ity, we limit the number of curves shown in figure [}] and 
omit the results for the ACDM model. For this model the 
scatter is similar to the OCDM model, with a somewhat 
higher amplitude on large scales and a somewhat lower 
amplitude on arcminute scales. The scatter in the value of 
te increases for larger aperture sizes. The value of te is 
independent of aperture size, and as a result the relative 
contribution of large scale structure to the error budget 
increases with increasing aperture size. Also the scatter is 
larger for fainter limiting magnitudes, i.e. higher source 
redshifts. 

Figure ||b shows the scatter for the (^-statistic. Because 
both (Clss) 1 ^ 2 an d 0) decrease with increasing aper- 



ture size the interpretation of figure is not as straight- 
forward as is the case for figure ||a. The effect of aperture 
size on the accuracy of weak lensing mass estimates is 
discussed in detail in section 4.3. 

Comparison of the results for the SCDM and OCDM 
cosmologies indicates that the differences are fairly small. 
We found that the ACDM cosmology gives similar esti- 
mates. Although we consider these different cosmologies 
throughout the paper (and plot the results for the OCDM 
and SCDM models) we note that the results hardly de- 
pend on the assumed cosmology. 

4. Results 

We examine the uncertainty in the weak lensing mass es- 
timate caused by uncorrelated large scale structure. The 
results are derived for a cluster with a SIS mass profile, 
and a velocity dispersion of a = 1000 km/s, typical for 
rich clusters that have been studied in the past. We study 
the dependence on cluster redshift, where we discuss the 
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results for low redshift clusters in more detail. We also 
examine how the accuracy of the mass measurements de- 
pends on aperture size. 

4-1- Dependence on cluster redshift 

First we examine how the mass estimates for clusters at 
different redshifts are affected when the lensing signal is 
measured out to a radius of 1.5 Mpc. 

For the SIS model we use the measurements in an an- 
nulus from 0.225 to 1.5 Mpc (R = 1.5 Mpc, 
a = 0.15). In the case of the (—statistic we use an in- 
ner radius of 9\ = 1 h$Q Mpc and an outer radius of 
02 = 1.5 h^Q Mpc (a = §). The value of the (-statistic 
only provides a lower limit to the cluster mass, and one 
way to circumvent this is to estimate the average surface 
density in the annulus around the aperture, for instance 
by using the results from a model fit. However, this model 
fit is also affected by large scale structure. In general the 
region where the model is fitted to the data largely coin- 
cides with the annulus used to compute the value of the 
(■-statistic. Thus if the value of the (-statistic is increased 
by additional structures, the mass from the model will also 
be higher. Although the uncertainty in the model param- 
eters and the (-statistic are not completely correlated, we 
will assume they are. This gives a fair upper limit on the 
uncertainty in the estimate of R 



LSS LSS 



r LSS 



(1 + a)6 2 



(25) 



The results are presented in figure Panels c and f 
show respectively the value of the Einstein radius and R 
as a function of cluster redshift. We show the results for 
the SCDM (thick lines) and the OCDM (thin lines) cos- 
mologies. The solid lines are for a limiting magnitude of 
R — 24 and the dashed lines for R = 26. Panels b and e 
show the scatter in the aperture mass statistics introduced 
by the large scale structure. Finally the ratios of the scat- 
ter and the lensing signal are presented in figure a and d. 

For clusters at low redshifts (z < 0.2) the relative im- 
portance of large scale structure for the ACDM model is 
similar to that of the SCDM model. For clusters at high 
redshift, or deep observations, the results for the ACDM 
are close to the OCDM results. 

Comparison of figures [|a and d shows that the relative 
importance of the distant large scale structure is smaller in 
the case of the SIS model fit. For this mass estimator the 
uncertainty in the cluster mass estimate is about 6% for 
rich clusters at intermediate redshifts. The uncertainty in 
the mass is large for high redshift clusters when the limit- 
ing magnitude of the sources is rather bright (solid lines) . 
However, in practice, to detect a significant lensing signal 
of a high redshift cluster deep observations are needed to 
ensure a sufficiently high number density of sources. As a 
result the importance of large scale structure is minimized 




lim 



Fig. 5. (a) Contour plot of the signal-to-noise ratio for a 
SIS model fit to the tangential shear of a cluster of a = 
1000 km/s as a function of limiting magnitude Ru m and 
cluster redshift z c . The noise includes the contribution of 
the intrinsic shapes of the sources and the contribution by 
large large scale structure for the OCDM cosmology, (b) 
Contour plot of the signal-to-noise ratio when the noise 
caused by large scale structure is ignored. The interval 
between adjacent contours is 2. 



(dashed lines). In the case of nearby clusters, the presence 
of structure along the line of sight results in a relatively 
large uncertainty in the cluster mass (see section 4.2). 

The uncertainty introduced by large scale structure is 
much larger when the (-statistic is used. This is not sur- 
prising as one only uses the shear measurements at rela- 
tively large radii, where the cluster mass surface density 
has become fairly low. In addition to this annulus, the 
SIS model fit uses data at much smaller radii, where the 
cluster signal is much higher. 

The noise in weak lensing analyses is generally domi- 
nated by the intrinsic shapes of the source galaxies, and 
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not by the noise in the large scale structure. Both sources 
of error depend differently on the limiting magnitude of 
the background galaxies, and as a result the total noise is 
a rather complex function of cluster redshift and limiting 
magnitude, as is demonstrated in figure for the OCDM 
cosmology. 

Figure |^b shows the the signal-to-noise ratio for a SIS 
model fit to the tangential shear of a cluster of a — 1000 
km/s, when the contribution of large scale structure is ig- 
nored. It should be noted here that we assume that the 
shapes of the sources have been measured accurately. In 
reality the shapes of very faint galaxies are difficult to de- 
termine. From figure ^|b one would expect to be able to 
measure accurate masses for low redshift clusters. How- 
ever, the results change significantly when we include the 
contribution from large scale structure as is shown in fig- 
ure [|a. The change in signal-to-noise ratio is small for 
high redshift clusters, where the shapes are still the ma- 
jor contributor to the noise. The true signal-to-noise ratio 
is lowered significantly in the case of deep observation of 
clusters at intermediate redshifts, but the largest differ- 
ence is seen for clusters at low redshifts. 

The signal-to-noise ratios for the ACDM model are in- 
termediate to those found for the OCDM and SCDM mod- 
els, resembling the OCDM better for high redshift clusters. 
However, at high redshifts the differences between the var- 
ious models become less important (they are fairly small 
anyway), because the noise is dominated by the intrinsic 
shapes of the sources. 

4-2. Low redshift clusters 

Figure ^|a demonstrates that large scale structure intro- 
duces a large uncertainty in the mass measurement of low 
redshift clusters. To date not many low redshift clusters 
have been studied through weak lensing. This is because 
the amplitude of the lensing signal peaks at intermedi- 
ate redshifts, and the signal is extremely small for nearby 
clusters. Furthermore their angular sizes were too large for 
the CCD cameras. Panoramic CCD cameras have become 
available recently, and a number of low redshift clusters 
have been imaged for weak lensing purposes. For instance 
Joffre et al. (2000) have analysed the mass distribution of 
Abell 3667 (z = 0.055). 

Although the lensing signal of low redshift clusters is 
small, there are several reasons to study these systems. 
First of all, they can, or have been studied extensively 
by other techniques. Also, the uncertainty in the redshift 
distribution of the sources is almost non-existent, and one 
can use relatively large background galaxies. 

In this section we examine how the weak lensing mass 
measurement of the Coma cluster (z — 0.023) is affected 
by structures at higher redshifts. A weak lensing analysis 
of the most nearby rich cluster of galaxies is challenging, 
but could provide important insights in the cluster dy- 
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Fig. 6. (a) Sources of noise in the measurement of the 
Einstein radius te from a SIS model fit. The figure shows 
the noise caused by the large scale structure (dashed lines) 
and the intrinsic shapes of the sources (dotted lines) as a 
function of limiting magnitude R\i m of the sources. The 
thick lines correspond to the SCDM model and the thin 
lines to the OCDM cosmology. The results for the ACDM 
model are similar to the SCDM model, (b) Expected 
signal-to-noise ratio as a function of limiting magnitude 
of the sources for the Coma cluster. For limiting magni- 
tudes R > 24 the signal-to-noise ratio of the measurement 
of the Einstein radius no longer increases, because of the 
noise from the large scale structure. The results for the 
ACDM model lie between the drawn curves. 



namics. We assume a velocity dispersion of 1000 km/s for 
the cluster (e.g., Colles & Dunn 1996). 

In figure |]a we plot the sources of noise in the measure- 
ment of the Einstein radius as a function of limiting mag- 
nitude of the sources. The figure shows the contribution 
from the large scale structure (dashed lines) and the intrin- 
sic shapes of the sources (dotted lines). For the intrinsic 
shapes of the sources we have assumed (7^) 1 ^ 2 = 0.25, 
similar to what is observed in deep HST images (e.g., 
Hoekstra et al. 2000) Reaching the depth of the HDFs one 
would expect to obtain a signal-to-noise ratio of more than 
20. However, figure ^a demonstrates that the total uncer- 
tainty in the measurement of the Einstein radius reaches 
a minimum value for a limiting R = 24. 

Figure ^|b shows the expected signal-to-noise ratio as 
a function of limiting magnitude, and it demonstrates 
that observations with Ry nll > 24.5 no longer improve the 
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Fig. 7. As figure |6| but now for the C~ statistic. The noise 
introduced by the large scale structure results in a con- 
stant signal-to-noise ratio for limiting magnitudes R > 25. 
Note that the achieved signal-to-noise ratio for the (- 
statistic is lower in comparison to the SIS model fit. 



signal-to-noise ratio of the mass measurement. The max- 
imum signal-to-noise ratio that can be reached by fitting 
a SIS model to the observed tangential shear around the 
Coma cluster is ~ 5 for the SCDM model, and ~ 7 for 
the OCDM model. The results for the ACDM model lie 
between the curves for the other two models. 

In figure [?] we show how the accuracy of the mass mea- 
surement is affected by large scale structure when the C - 
statistic is used. Comparison with figure ^ shows that the 
SIS model fit yields a better signal-to-noise ratio. 

4-3. Large apertures 

The intrinsic shapes of the sources result in a statistical 
uncertainty in the measurement of the Einstein radius. If 
the SIS model is fitted to the measurements in an annulus 
from Rq out to R, the error in te is given by 



nn ]n(R/Ro) 



Cgal, 



(26) 



where n is the number density of sources, and er ga i is the 
uncertainty in the measurement of the tangential distor- 
tion for a single galaxy. 

From equation ^ we find that this contribution to the 
error budget decreases slowly for increasing aperture sizes. 



Fig. 8. (a) Error on the measurement of the Einstein ra- 
dius te as a function of aperture size, because of the 
large scale structure (dashed lines) and the instrinsic 
shapes of the sources (dotted lines) . The SIS model is fit- 
ted to the tangential shear in an annulus ranging from 
Ro = 0.225 Mpc out to R out . The thick lines give the 
results for the SCDM model, whereas the thin lines corre- 
spond to the OCDM model, (b) The expected signal-to- 
noise ratio as a function of aperture size. The curves for 
the ACDM model are comparable to the results for the 
OCDM model. The signal-to-noise ratio does not improve 
for aperture sizes i? ou t > 1 h$Q Mpc. 



However, figure shows that the noise contribution from 
large scale structure increases rather rapidly with aperture 
size. 

Figure |^a shows how the contributions of the intrin- 
sic shapes of the sources, and the large scale structure 
depend on the aperture size. The results are for a clus- 
ter at a redshift z — 0.5, and the lensing signal has been 
measured from R = 0.225 Mpc out to a radius i? ou t- 
The uncertainty in te caused by the intrinsic shapes of the 
sources decreases slowly with aperture size (dotted curve) , 
whereas the contribution by the distant large scale struc- 
ture increases (dashed line). Figure ||b shows the resulting 
signal-to-noise ratio. The signal-to-noise ratio is virtually 
constant for outer radii larger than ~ 1 h§ Mpc. Measur- 
ing the lensing signal in larger apertures does not improve 
the accuracy of the cluster mass. One has also to keep in 
mind that we have assumed that the mass profile is well 
described by a SIS model, which might not be true in the 
outer parts of clusters. 
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For clusters at higher redshifts the curves look similar, 
but the maximum signal-to-noise ratio is reached at some- 
what larger (physical) aperture sizes. The noise caused by 
the large scale structure dominates the error for low red- 
shift clusters, and for these systems, the signal-to-noise 
ratio decreases slowly with increasing aperture size. 

So far we have used an inner radius for the SIS model 
fit of i?o = 0.225 h^Q Mpc, but one cannot always fit 
the model at such small distances from the centre. For 
instance in the case of the z — 0.83 cluster MS 1054-03 
Hoekstra et al. (2000) showed that the mass distribution 
in the centre is rather complex and the substructure lowers 
the azimuthally averaged tangential shear. Consequently 
they fit the SIS model from - 0.7 out to - 2 Mpc, 
which corresponds to a = 0.35. 

We also have examined how the accuracy of the mass 
estimate changes when the SIS model is fitted to the 
data within a given aperture, while the inner radius is 
increased. In this situation both the noise caused by the 
shapes of the sources and the contribution by the large 
scale structure increase. Thus small values for the inner 
radius improve the signal-to-noise of the mass measure- 
ment, but one should keep in mind that substructure in 
the cluster centre complicates a safe choice for the inner 
radius. 

For the (-statistic the noise contribution caused by the 
shapes of the galaxies is given by 



1 



tt(1 - a 2 )9f 



-CTgal- 



(27) 



However, we want to estimate R(< 1 h% c fMpc), and 
thus we have to add the average surface density in the 
annulus from 1 to R ou t Mpc. As before we estimate 
R(l < r < Rout) from a SIS model fit to the tangential 
shear. The errors in the (— statistic and the Einstein radius 
from the model fit are combined according to equation 25, 
because the errors are correlated. We note that the size 
of the correction depends on the assumed mass profile, 
but for increasing values of i? ou t the correction becomes 
smaller (i.e., R(l < r < R out ) becomes smaller). 

The resulting noise estimates are presented in fig- 
ure 0a. The corresponding signal- to- noise ratio of R(< 
1 h$ Mpc) for a cluster at a redshift z = 0.5 is presented 
in figure |9(b. It shows that using a larger i? ou t improves 
the accuracy in the measurement of R(< 1 h§ Mpc). 

A similar calculation shows that the contribution of 
the large scale structure becomes more important when 
one tries to measure the average surface density in larger 
apertures. For instance the uncertainty in k(< 2 h^ Mpc) 
almost doubles compared to R(< 1 /i^Mpc). 

Attempts to measure the average cluster mass profile 
out to large projected distances from the cluster centre 
will also be affected by large scale structure, given the 
fact that it is already difficult to measure the total mass. 
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Fig. 9. As figure ra, but now we show the measurement 
error on R(< 1 h§ Mpc) as a function of Rout when the 
(■—statistic is used to estimate the mass. The average sur- 



face density in the annulus from 1 to R Q 



'50 



Mpc is 



estimated from a SIS model fit to the observations. The 
thick lines give the results for the SCDM model, whereas 
the thin lines correspond to the OCDM model, (b) The 
expected signal-to-noise ratio of R(< 1 /i^Mpc) as a func- 
tion of aperture size. The signal-to-noise ratio improves 
with increasing i? out . If the cluster mass profile is well 
described by a SIS model, the model fit gives the better 
results. 



Unfortunately it will be difficult to quantify the effect of 
distant large scale structure on such measurements. 

5. Conclusions 

The observed weak lensing signal is sensitive to all mat- 
ter along the line of sight, and as a result this additional 
matter acts as a source of uncertainty in mass estimates of 
clusters of galaxies. In this paper we have investigated the 
effect of distant large scale structure, uncorrelated with 
the cluster mass distribution. It can be treated as an ad- 
ditional source of noise. 

Given a cosmological model it is possible to compute 
the statistical uncertainty in the weak lensing mass mea- 
surement caused by the distant large scale structure, and 
quantify its contribution to the total error budget. We 
have considered two methods to estimate the mass from 
a weak lensing analysis. One way is to fit a SIS model to 
the observed lensing signal. The other method uses the (- 
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statistic (e.g., Fahlman 1994) to estimate the mass. Both 
methods give similar results. 

We find that the importance of matter along the line 
sight is fairly small for rich clusters (a = 1000 km/s) at in- 
termediate redshifts, provided that the bulk of the sources 
are at high redshift compared to the cluster. If the lens- 
ing signal is measured out to 1.5 h§ Mpc from the clus- 
ter centre, the typical la relative uncertainty in the mass 
cause by the large scale structure is about 6%. 

In certain situations the effect of large scale structure 
is more important. For nearby clusters, such as the Coma 
cluster, background structures introduce a considerable 
uncertainty in the mass. Considering only the noise caused 
by the shapes of the sources used in the weak lensing anal- 
ysis one expects to achieve very accurate weak lensing 
mass measurements of nearby clusters, but we find that 
the large scale structure limits the maximum achievable 
signal-to-noise ratio to ~ 7, irrespective of the depth of 
the observations. 

We also have examined how the accuracy of the weak 
lensing mass estimate changes when the lensing signal is 
measured out to large projected distances from the cluster 
centre. In the case of the SIS model fit we find that the 
signal-to-noise ratio does not increase significantly when 
one includes measurements at radii larger than ~ 1 h^ 
Mpc. When the (^-statistic is used to estimate the mass 
within a 1 h§ Mpc radius aperture the signal-to-noise 
ratio does increase with i? ut- If the cluster mass profile 
well described by a SIS model, the model fit results in a 
better signal-to-noise ratio. 

The effect of large scale structure is important in the 
outskirts of clusters, where the cluster signal itself is low. 
In this paper we have quantified the uncertainty in the 
measurement of the cluster mass, but distant large scale 
structure will also introduce noise in studies that try to 
constrain the mass profiles of clusters at large projected 
distances. Furthermore studies of the 3-D structure of 
clusters by combining weak lensing, X-ray, and S-Z data 
(e.g., Zaroubi et al. 1998) will be affected by uncorrelated 
large scale structure. 

Photometric redshifts of the sources can be used to de- 
crease the effect of large scale structure, either by model- 
ing the contribution of large scale structure to the lensing 
signal, or by selecting sources in an optimal redshift range. 
This can be only done at the expense of a large amount 
of additional data. Because the effect of large scale struc- 
ture can be treated as an additional (statistical) source of 
noise. Thus an alternative procedure is to average the re- 
sults of an ensemble of clusters, thus decreasing the noise 
on the (average) cluster mass estimate. 
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